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Abstract 

■ The Schrodinger equation for an electron near an azimuthally symmetric curved surface E in the 

presence of an arbitrary uniform magnetic field B is developed. A thin layer quantization procedure 
is implemented to bring the electron onto E, leading to the well known geometric potential Vc oc 
h? — k and a second potential that couples An, the component of A normal to E to mean surface 
. curvature, as well as a term dependent on the normal derivative of An evaluated on E. Numerical 

results in the form of ground state energies as a function of the applied field in several orientations 
are presented for a toroidal model. 
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1. INTRODUCTION 



Nanostructures with novel geometries have become the subject of a large body of exper- 
imental and theoretical work |j H S i, 1 ©, 0, 1, S El ED, El El El, El El E3 El 
EH, HI EH HI 123]. Many of the fabricated structures exhibit curvature on the nanoscale 
making once purely theoretical investigations of quantum mechanics on curved and reduced 
dimensionality surfaces relevant to device modelling. Because magnetic field effects prove 
important to Aharanov-Bohm and transport phenomena, the inter play of an applied field 
with the curved regions of a nanostructure 0, H H M, El, El l20l 124 l25j (particularly if 
the structure has holes) may be critical to a complete understanding of the physics of a 
nanodevice element. 

In a previous work || the Schrodinger equation for an electron on a toroidal surface T 2 
in an arbitrary static magnetic field was developed and numerical results obtained. There 
however, the electron was restricted ab-initio to motion on T 2 , which pre cluded the appear- 
ance of the well known geometric potential Vc oc (h 2 — k) , [H |27], l28l 1291 f30L l3ll 1321 . I33T l34l 
I35L l36l . H3, IH, HI HI Ell EH with h, k the mean and Gaussian curvatures, that arises via 
thin-layer quantization |43|. Furthermore, if the degree of freedom normal to E (labelled q 
in everything to follow) is included, then the component of the vector potential An normal 
to E couples to the normal part of the gradient. Given that Vc is generated from differenti- 
ations in q and the requirement of conservation of the norm, it is not surprising that other 
curvature effects would follow from inclusion of A. This work is concerned with determining 
the effective potential that arises from the A^d/dq operator, and through numerical calcu- 
lation with a simple but realistic model, gauging its influence on the single particle spectra 
and wave functions as a function of field strength and orientation. 

The remainder of this paper is organized as follows: Section 2 develops the Schrodinger 
equation 

1 f h \ 2 

— (-V + eA) ty = E^ (1) 

2m \ i J 

for an electron allowed to move in the neighborhood E. In a preprint version of this work 44] 
differential forms were employed to represent Eq. (1), but here more conventional language 
is used. Section 3 briefly reviews the procedure by which Vc is derived and uses a similar 
but not identical methodology to reduce the A^d/dq operator to an effective potential 
expressed entirely in surface variables. Section 4 employs the formalism presented in section 
3 to calculate spectra and wave functions as a function of field strength and orientation for 
a toroidal structure, and section 5 is reserved for conclusions. 



2. CURVED SURFACE SCHRODINGER EQUATION 

The development of the Schrodinger equation in three dimensions on an arbitrary mani- 
fold generally yields a cumbersome expression. To remove some complexity in what ensues, 
an azimuthally symmetric surface E with q the coordinate that gives the distance from E 
will be adopted. 

Let e p , e^, e z be cylindrical coordinate system unit vectors. Parameterize E(p, 0) by the 
Monge form 

r(p, (f)) = p e p + S(p) e z (2) 

with S(p) the shape of the surface. Points near the surface S(p) are then described by 

x(p,0,g) = r(p,<j>) + q e n (p,0) (3) 

with e„ everywhere normal to the surface and to be defined momentarily below. The differ- 
ential line element of Eq. (3) is 

(ix = dr + dq e n + q de n . (4) 
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After some manipulation along with a list of expressions to be defined below, Eq. (4) can 
be rewritten as (the subscript on S denotes differentiation with respect to p) 

(5) 



dx = Z(l + k\q)eidp + p(l + k 2 q)e ( f ) d(J) + dqe Ti 
= Zfieidp + pfa%6&<\> + dqe n 



with 



1 + s 2 



ei 



— (e p + S p e Z/ 



and the principle curvatures 



S, 



pp 



ko = — 



Z 3 
pZ' 



The metric elements can be read off of 

dx 2 = Z 2 f 2 dp 2 + p 2 f 2 d(j) 2 + dq 2 



(6) 
(7) 
(8) 

(9) 
(10) 

(11) 



from which the Schrodinger equation can be determined, but since the minimal prescription 
will be employed as per Eq. (1) it proves convenient to use the gradient 



1 d 



}\Z dp f 2 p 
instead. Eq. (1) can be rearranged to 



I d_ d_ 



(12) 



2 



+ 



Id lid 2 



2lz*fldp* hhZpdp hf 2P 



o 2 



ki ki\ d 
fi fi J dq dq 2 



pk\k 2 h gk lp ^ d_ 



ft 



2Xi 



A l d A, d 



fi ^ ft J dp 

X 2 (A 2 1 + Al + A 2 N )+ h2 



+ 



fxZdp f 2 p< 
2Em~ 



a 9 
dq 



* = 



(13) 



with A = e/h and Aj = A • ej. 

While Eq. (13) describes the general case for an azimuthally symmetric geometry, it can 
be simplified substantially when considering the q — ► limit, or as dubbed by Golovnev [43], 
with thin layer quantization. The procedure entails first performing all q differentiations in 
accordance with the intuitive notion that the kinetic energy in a thin layer is large, then 
setting q = everywhere. Eq. (13), leaving the ^-differentiations intact and setting q = 
everywhere save the A^ term, cleans up to 



Id 2 lid Id 2 d d 2 
Z 2 dp 2 Z pdp p 2 d(f) 2 dq dq 2 
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d_ + A t 
, dp P 



g=0 



— + A — 

d(j) N dq 



+ 



q=0 



-X\A\ + Al + Al) 

with the mean curvature h above given by 

h = ^(k 1 + k 2 ) 

The Gaussian curvature which will appear later is 

k = fci&2- 



2Em 



^ = 



(14) 



(15) 



(16) 



3. DERIVATION OF GEOMETRIC POTENTIALS 



The geometric potential Vc(p) is found by reducing Eq. (13) further by the well known 
procedure of assuming a suitable confining potential in the normal direction V n (q) and 
demanding conservation of the norm in the q — ► limit. The latter requirement is generally 
imposed via assu ming a separation of variables in the surface and normal parts of the total 
wave function |2f| UA EE HO, HO, EO, EH 

*{p,M ^Xs{p,<P)XN(q) (17) 



and imposing conservation of the norm through 



\y(p,<P,q)\\l + 2qh + q 2 k)di;dq ^ \xs(pA)\ 2 \XN(q)\ 2 dIldq (18) 

or equivalently 

* = XsXtv(1 + 2g/i + g^)- 1 / 2 = XsXnG~ 1/2 . (19) 

Inserting the rightmost term in Eq. (19) into Eq. (13) and subsequently taking q — ► 
reduces the g differentiations there to 

d d 2 d 2 2 
"^<9g <9g 2 <9g 2 ^ ^ 



Separability of the surface and normal variables in the non-magnetic part of the Hamiltonian 
is manifest. 

Since An can be a function of q while h is not, it is not immediately apparent that 
separability of the Schrodinger equation in the surface and normal variables is preserved 
by the AN(p,4>,q)d/dq operator. Rather than apply the identical procedure above to this 
term, instead integrate out any g-dependence with some reasonable ansatz for Xn{q), say 

a normalized hard wall form Xn(c[) = y/^/L simrq/L. While equivalent to the method 
summarized above in Eqs. (17)-(20), it assists in establishing conditions on A^ that will 
prove useful later, and previous work [45] has demonstrated that to a good approximation 
this procedure is justified even if a particle is not strongly confined to a region near E. 
Write 



d Xn(q) 
dq G-V2 



Gdq 



(21) 
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which is equivalent to 



XnWQ 1 + qk)G 1 A N (p,(p,q)dq+ / XN(q)A N (p, (p, q)x' N {q)dq. (22) 



I = - 



Consider the left hand integral II, and assume (suppressing surface arguments) A N (q) ~ 
a + aiq + expanding out the arguments and noting that each power of q when integrated 
picks up a power of L that will vanish as L — > 0, the result in this limit is 



1 L — 



-hA N (p,(f),0) / x 2 N(q) d( i 



(23) 



i.e., the effective potential arising from I L is proportional to — hA N (p, 0, 0). 

Turning to the second integral Ir, write Xn(q)x'n(q) as 9/dq[x%(q)/2] and perform an 
integration by parts. The surface term vanishes so that 



L 2 / s9A N (p,(f),q) 
XnW Q-q q - 



(24) 



Again taking A N (q) ~ a + a\q + ... allows Eq. (24) to be expressed in the q — > limit as 

dq. (25) 



ru ~ 
2 R — 



L 2 / ^9A N (p,(f),q) 
XnW- 



dq 



q=0 



Eqs. (23) and (25) combine to give an effective surface potential (with constants appended) 

(26) 



h(p)A N {p, 0, 0) + — 



g=0J 



There are two points that should be addressed with further explanation. First, assuming a 
nonsingular series expansion for A^ is reasonable since A^{p, 0, 0) is the value of the vector 
potential on the surface; it should certainly be physically well behaved there but need not 
vanish. Secondly, while a hard wall form has been assumed for the trial wave function (a 
Gaussian works as well), it is a reasonable conjecture that the arguments made above are 
independent of the choice of XN{q) as long as there is negligible mixing amongst states in 
the q degree of freedom. 



4. NUMERICS 



The formalism developed above is now applied to calculate the spectrum and eigenfunc- 
tions for a toroidal structure in a uniform magnetic field of arbitrary orientation inclusive 
of geometric potentials. 

A convenient choice to parameterize points near a toroidal surface T 2 of major radius R 
and minor radius a is [46] 

x(0, 0) = W(6)e p + a sin6»e z + qe n (27) 

with 

W = R + a cos9, (28) 

and 

e n = cos8e p + sin6>e 2 . (29) 
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The differential line element is 

rfx = a(l + kjq)e e d9 + W(9)(l + kjq)e (j) d(f) + e n dq (30) 

= a q (q)e e d9 + W q (Q, q)e (j) d<j) + e n dq 
with e e = — sin#e p + cos#e z and the toroidal principle curvatures 

kl = (31) 



*-m (32) 

The gradient that follows from Eq. (30) is 

^ 1 d 1 d d 

a q (q)d9 v W q (9,q)d<j) dq 

The symmetry of the torus allows an arbitrary static magnetic field to be taken as 

B = Bii + B k. (34) 

In the Coulomb gauge the vector potential A(9, 0) = |B x r expressed in the geometry of 
Eq. (27) is 

A(9, (f), q) — ^ [i?isin0(i?cos^ + a g )eg + (-BoW 7 ? — -BiOq sin#cos0)e^ 

+i? 1 _Rsin0sin 9e n ]. (35) 

In this case dA^/dq = so only the first term in Eq. (26) will contribute to V^ ag (9 , <j)) . 
The Schrodinger equation in the q — > limit inclusive of geometric potentials can be written 
(spin will be neglected) in a compact form by first defining 

a = a/R 

F(9) = 1 + a cos9 

70 = B 7iR 2 

71 = .BxTri? 2 

7w = — 

e 

7o 

To = 

In 

In 
2mEa 2 

e = - 



h 2 



after which Eq. (1) may be written 



d 2 a sin 9 d a 2 d 2 1 iar x . . (1 + 2a cos9) 

W ~ ~F~Wd9 + F\9) <90 2 + 4F^9) + — sm9sm( P pffi 
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T a 2 - — — sin#cos</> — + iar^uKpia + cos#) — 
F(8) J d(p ' oO 

r 2 a 2 F 2 (6) T 2 a 2 F 2 (6) ( . 2j a 2 sin 2 0\ r O 7ia 3 F(0) . „ / 

4 - ( sm + -F^rJ + — 2"^^«^_ 



* = (36) 



=5> if T * = (37) 

with the fourth and fifth terms of Eq. (36) being proportional to Vc and V™ ag . 

To obtain solutions of Eq.(36) a basis set expansion orthogonal over the integration 
measure dJ(9) = F(8)d8d(p may be employed for a given a. Here R will be set to 500 A in 
accordance with fabricated structures [T3, |23|, H3, |48[ (for an _R = 500A torus r = .263-Bo)- 
With a = 250A, a = 1/2, a value which serves as a compromise between smaller a where 
the solutions tend towards simple trigonometric functions and larger a which are less likely 
to be physically realistic. 

There are two options available for basis functions, one being 



Xr, 

and the other being 




Xt(0,4>)= Jn \ e** (39) 



with f n (8),g n (9) even/odd orthonormal functions labelled by +/- respectively, constructed 
by a Gram-Schmidt (GS) procedure |4j| over dJ(8) using cos n6, sin n6 as primitives. Here 
the latter approach will be adopted in order to facilitate comparisons to a related work 0] 
and because of the ease by which the Hamiltonian matrix elements 

Ununv ~ {Xnu\-^r\Xnu) (40) 

can be evaluated analytically. The basis comprises six GS functions of each ^-parity and 
five azimuthal functions spanning —2 < v < 2 per ^-function for a total of 60 basis states. 
The resulting 60 x 60 Hamiltonian matrix blocks schematically into 




from which eigenvalues and eigenfunctions are determined. Since the concern here is with 
the ground state eigenfunctions only a few GS states prove relevant; they are fo(9) = .3987, 
f^O) = .6031 cosfl - .1508 and g x (ff) = .5642 sinfl. 

Figs. (1-3) plot the ground state energy e as a function of the magnitude of flux r for 
three field orientations. Each figure displays results for Eq(t) with Vc and/or V mag switched 
on or off. In Fig. 1 the field is oriented along the z-axis. In this case there there is no 
component of A normal to the surface so only two curves are evidenced. The curves are 
qualitatively similar and the effect of Vc is to smooth out the Vc = curve at r ~ 1 and 
shift it downward by an overall constant. It should be noted that although the curves are 
similar, persistent current effects depend of the smoothness and shape of the of e{r) curves 
so that even fine details can prove important. Table I shows the evolution of the ground 
state wave function Xg(9,4>) for several r; because the field is oriented along the z-axis, r 
also measures the flux through the toroidal plane. Both Vc and V ma9 are zero so that the 
evolution of Xg{9, 4>) is due only to changes in field strength. 
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Fig. 2 shows results for a field orientation oriented tilted | radians relative to the toroidal 
plane with the magnitude of r plotted on the horizontal axis. The divergence of the two 
lower curves illustrates the influence of V mag on the spectrum. The curve inclusive of V mag 
begins to involve excited states in both the 9 and azimuthal degrees of freedom (see table 
II). While the admixture of excited states would tend to raise Eq, the interaction is strong 
enough to pull the e (t) down as the applied field increases. 

In Fig. 3 the magnetic field is situated parallel to the toroidal plane along the x-axis. 
As would be anticipated, there is no structure in the Vc = V mag = curve since no flux 
penetrates the plane of the torus. This trend obtains also when Vc, V mag ^ 0. The effect 
of V mag becomes substantial very quickly both in Eo(t) and on Xg(#>0) (table III). In this 
case there is not even qualitative agreement between results with Vc ^ with V mag omitted 
compared to when it is included. 



5. CONCLUSIONS 

This work presents a method to reduce the A^d/dq term appearing in the Schrodinger 
equation for an electron near a two-dimensional surface in an arbitrary static magnetic field 
to a geometric potential V ma9 written entirely in terms of surface variables. This potential 
can appreciably modify energy vs. magnetic flux curves as well as surface wave functions 
considered here. 

In the context of real structures, the practical utility of deriving geometric potentials 
lies in reducing three dimensional problems to two-dimensional ones. Earlier work 45] on a 
simple nanoscale model has shown that solutions of an ab-initio two dimensional Schrodinger 
equation do not adequately approximate thin layer three-dimensional solutions on a curved 
space unless Vc is included in the Schrodinger equation. From this perspective, geometric 
potentials of the form discussed here should be considered effective potentials that must be 
included in the modelling of curved structures in order to achieve a complete description of 
the object. 
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Table I. Ground state wave functions for B = B h at integer values of r. The arguments in square 
brackets indicate if Vc/ V mag are switched off or on. Only the dominant terms are shown. 





r = 


r = 1 


r = 2 


XG [off, off] 
Xc[on,off] 
Xc[on,on] 


1 

-.968/o - -244A 
-.968/q - .244A 


~ 1 

-.957/q + .254/x 
-.957/q + -254/i 


(-.969/o + .245/i)e-^ 
(.987/o - .158/i)e-^ 
(.987/q - .158/i)e-^ 



Table II. Ground state wave functions for B = Bq(\ + k)/\/2 as per table I. For this field configu- 
ration An(9, (p) / 0. 





r = 





r = 


-- 1 


T = 2 


XG[0ff,0ff] 


1 




-.989/o - 


.115flFie -< * 


(-.928/o + .159/i)e-^- 












51 (.287- .145e~^) 


X G [on,off] 


.968/q - 


.244/! 


.961/o - 


- .252/i 


.932/o - .270/i 


Xc[on,on] 


•968/q - 


.244/i 


(.957/o - 


•232/i)+ 


/o(.909e"^ - .126)- 








#i(.094e^ - 


- .me-**) 


ff i(.173e-^-.351) 



Table III. Ground state wave functions for B = B\i at integer values of r as per table I. An(9, (j)) 
is nonzero, and there is no flux through the toroidal plane. 





r = 


r = 1 


r = 2 


XG [off, off] 
X G [on,off] 
Xc[on,on] 


.968/q - -244/1 
.968/q - .244/i 
.968/o - .244/i 


.978/o + .279%sin0 
.964/o - .218/i 
-.954/o + .178/i + .320%sin</> 


.894/o + .133/i - .552%sin0 
.941/o - .132/i + .4O3w/isin0 
.869/o - .250/icos</> - .314% sin(/) 



11 




FIG 1. £o as a function of r for B = Bok. Diamonds correspond to Vc = V mag = 0, stars to 
V c / 0, V ma 9 = and squares to V c + 0, V ma 9 / 0. 
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FIG 2. £o(r) for B = B (i + k)/V2. Diamonds correspond to V c = V ma9 = 0, stars to V c ± 0, 

yrnag = q and gquares t O Vb ^ 0, V™ 9 / 0. 
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FIG 3. £ (r) for B = B^i. Diamonds correspond to V c = V ma9 = 0, stars to V c + 0, V ma 9 = 
and squares to V c + 0, y mafl / 0. 
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